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I. INTRODUCTION 



Variational principles, well established in nonrelativistic quantum problems, develop 
their true power by setting a measure for the best approximate solution within some 
parametrization of a trial space. This space is made up of wave functions ip in non- 
relativistic quantum mechanics, by statistical operators p in thermodynamics, and by 
actions S in the Feynman-Jensen formulation. The above object if) can always be un- 
derstood to be the ground state of some Hamiltonian H. Any statistical operator is also 
related uniquely to a Hermitean operator H (its "Hamiltonian"). Thus, in any of these 
three cases, we may talk about a theory to be parametrized and varied. The task is to 
find a class of theories coming reasonable close to the truth but to keep it simple enough 
for tractability. 

Each of the above three cases generalizes to field theory, if> becoming a wave functional, 
see e.g. flT], while H and S keep their meaning. For the formulations of the thermal 
variational principle, with H [2-10] or with S [5,10-15], we refer to the next Section as 
some part II of the Introduction. 

We focus on the application to pure gauge theory with particular interest in the hot 
gluon system. There are three fascinating aspects of this system. First, it distills out 
from the standard model by reducing the number of flavours to zero, while possibly still 
containing the whole non-Abelian mystery. Second, other than in the Higgs sector (the 
other distillate) its Lagrangian looks so simple : C = — F 2 /4 . Third, its high-temperature 
limit may be called understood, mainly based on the rearrangement [ jT5| , |T7| ] of diagrams 
at soft-scale outer momenta Q (key words: resummation, hard thermal loops). 

There are several basic problems and pitfalls at the very beginning when the varia- 
tional method contacts gauge fields. It is the subject of the present paper to make the 
calculus working at all. Hence no new results on the hot gluon system should be expected 
here. Instead, the known perturbative results are used to test our new variational con- 
struction. With this first step we hope to pave the road towards its presumedly powerful 
nonperturbative possibilities. 

Specifically on SU(iV) gauge fields under thermal variation, there is (to our knowl- 
edge) only the one (thus pioneering) paper of R. Manka in 1986. He studied pure 
non-Abelian gauge theory by using a free trial Lagrangian, namely the Abelian one for 
photons taken n-fold (n = N 2 — 1). The fields are identified, v4 non _Abeiian = -Atrial (at 
least in the high-temperature phase), and a constant transverse photon mass is taken as 
variational parameter. Note that this identification of pure, but interacting gauge fields 
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with free, but massive ones makes the trial theory nontrivial. With the longitudinal mass 
included as well as the 4-vertex (which both were neglected in 0), Manka conjectures 
that the perturbative results on masses, generated by the plasma, should be obtained 
from variation as well. Indeed, they should - but by this supposedly easy task we were 
led into all that follows. 

The following outlook reflects, to some extent, our individual path into the subject. 
Starting with the basic ideas of M just mentioned, we were more or less forced into one 
step (away from ||) after the other: 



1. The covariant gauge-fixing parameter a is reintroduced, and kept arbitrary, because 
all experience with, e.g., the damping puzzle of the gluon plasma tells us that a, if 
surviving in final results, is an ideal indicator for wrong physics. 

2. Both dynamical mass terms (transverse and longitudinal) are included as functions 
of momentum. This setup covers static screening as well as dynamically generated 
masses. The massive-photonic trial theory still keeps its Abelian gauge invariance. 

3. The functional integral formulation is applied. In passing, although our notation is 

Minkowskian (metrics H ) we actually always mean the Euclidean space. We 

only have to remember, at appropriate places, that the zeroth component A of the 
gauge field is i times a real field. 

4. The classical (or Feynman- Jensen) version of the variational functional is used, 
because it avoids difficulties in constructing the Hamiltonian to our higher-derivative 
trial-Lagrangian. As an intermediate result, the covariant gauge-fixing parameters 
of studied and trial theory become equal. 

5. The variational functional, if evaluated with the quadratic photonic trial theory as 
described, still depends on the (common) gauge-fixing parameter a (as also observed 
in [0). A way out is proposed by first rewriting the partition function of the theory 
studied such that its action becomes even in the gauge fields. This is called the "even 
version" in Sec. II C. 

6. As in the low-order perturbative treatments, and since we shall only reproduce its 
results, detailed renormalization is not (yet) required in this paper. Divergent terms 
can be separated from the finite thermal ones. Hence, the coupling g changes its 
meaning to be the running coupling in these thermal contributions. 

7. For a first application of the "even version", scalar electrodynamics is appreciated 



once more [19] to be an ideal toy model for the non- Abelian problem. The known 



2 



self-energies are put in by hand, but supplied with variable prefactors. Through 
variation, the latter become 1 indeed. 

8. In the non-Abelian case, the Faddeev-Popov determinant becomes part of the even- 
odd decomposition. The "even" functional works well, except for a (hopefully) minor 
detail at the end (concerning gauge-fixing dependence in higher order). 

9. For the explicit analysis just mentioned, the variational functional had to be ex- 
panded up to the third (partly fourth) g-power. This apparently inevitable recourse 
to ^-powers is a big disappointment. 

The paper is organized as follows. Section II on the formulations of the thermal 
variational functional is a continued Introduction. Especially the "even version" (the 
one that works) is introduced in Sec. II C. In Section III we follow the Feynman- Jensen 
version. It leads to unphysical results, but is, on the other hand, reasonable simple to 
introduce several technical details. Section IV treats scalar electrodynamics with the 
"even version" of the functional. In Section V on the gluon plasma, things start more 
involved but become very similar at the end. In Section VI the case of constant trial 
self-energies is discussed in terms of Debye screening and magnetic mass. Open questions 
are summarized in Section VII. Conclusions follow in Section VIII. Three Appendices 
cover details on the functional integral measures, on some normal integrals involved and 
on sum rules. 



II. THE THERMAL VARIATIONAL PRINCIPLE 

II A. Gibbs — Bogoljubov 

The extremal properties of thermodynamic potentials are known from text books 0, [| 
on statistical physics. In particular in the canonical ensemble (the only one considered 
in this paper), the free energy takes its minimum at equilibrium: F > F,. In its usual 
version, the thermal variational principle is identical with this modest inequality, if its 
left-hand side is detailed : 

V [H] = Tr (- [H.-H] J -Tln(Z) = min. . (1) 

The proof is given shortly. In ([[]) and in the following an index bullet refers to the system 
studied (at equilibrium), i.e. to the "hard problem" which one likes to learn about by the 
variational method. (3 — 1/T, Z — Tr ( e~^ H ) . Trial quantities carry no index, so H is the 
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element running through the trial space whose only restrictions are that (a) the spectrum 
of H is bounded from below and (b) H acts in the Hilbert space of H m . The formulation 
(HD is found e.g. as eq. (10.83) in or as eq. (20.37) in [|J. It is called Gibbs variational 
priciple in and Bogoljubov inequality in || |J. There is a natural application to the 
Heisenberg spin model, where mimimizing V yields the best Curie-Weiss Hamiltonian 
f|], thereby justifying the mean field procedure. 

For the proof of ([!]), we claim that one line suffices. It rests on the inequality — ln(x) > 
1 — x and on the irrelevance of operator-ordering under trace, Tr ln(AB) = Tr \n(BA). 
With any non-equilibrium statistical operator p, the line reads : 



F[p] = Tr (pH,) + T Tr (p In [ p ] ) = F, — TTr ( pin 



" 1 




-p. 


) 







> F. 



To the left, (0) starts with the non-equlibrium free energy in (E — TS)-form. The 
knowledge of p. at equilibrium, H, = — Tln(p.Z.), is used for the inner equality sign 
(F. = —Tin (Z,)). Finally, the right end has been simplified using Tr (p) = Tr (p.) = 1. 
In a whatsoever non-equilibrium state the system is, it has a statistical operator p with 
the three properties 1-trace, hermitecity and positivity. Thus (with the properties (a), 
(b) of H as stated above), its general form is p = e~^ H jZ. This makes (||) to become ([p), 
q. e. d. 

It is tempting to require that the trial theory be a solvable one (e.g. a free field 
theory). However, it must not. Imagine there was a small coupling e in the trial theory, 
and (for simplicity) only one variational parameter 77. Near its minimum, the functional 
would take the form V = a(e) +b(e) [i] — c(e)] 2 . Clearly, through perturbative expansion 
of V , the coefficients a, b, c as well as the position rj of the minimum would be obtained 
as power series in e. The parameter rj may be chosen to be the coupling e itself. 

We now turn to gauge field theory, governing a periodically repeated box of volume 
V and coupled to a thermal bath at rest with four-velocity U = (1,0). In the variational 
principle (|1|) H, and H are the Hamiltonians to a Lagrangian £. studied and a trial 
Lagrangian C, respectively. To count the same number of field degrees of freedom, one 



may either prepare the physical Hilbert spaces from the outset || || [15[ or work with 
extended spaces (corrected by ghosts). In the latter case the two gauge-fixings may be 
different. Adopting general covariant gauges, there is a gauge fixing parameter a, of the 
theory studied and an a of the trial theory. No final result is allowed to depend on either 
of them. 



The hot gluon system is described by the pure Yang-Mills Lagrangian 

1 1 / \2 

_ p a JP^ U a ___ 

4 ' 2a. 



C = -~F.^ a Fr a - ^- (d»A a Y + ic a d»Dfc b . (3) 
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with F,^ = d^A a u - d v Al + gf abc A\A c u and Df = 5 ab d^ - gf abc A c ^ . In the high- 
temperature limit the 1-loop contributions (hard thermal loops) are of relative order unity 
and therefore must be included in the true leading order Jl6], [L7| . The gluon propagator, 
resummed this way, may be written as 

r r fr y\ _ FT(Q) B"»(Q) D^(Q) (A] 

W) Qi-M t (Q) + Q*-M e (Q) +a > KV 

where R, B, D are members of the symmetric Lorentz matrix basis 

R-a B D Q- VoV r- Q° V + V °Q D _Q°Q ~ 
R-g-B-D , — , C " ^ Q 2 q ~ ~Q2~ (5) 

with V = Q 2 U - {UQ)Q = {-q 2 , -Q q ) • (6) 
The orthonormal properties of (0) are listed in [2~0|. Note that R and B are projectors: 

FHQ) Qu = , B""(Q) Q v = . (7) 
In (^), Mi = ZT t and = 77^ are the well known polarization functions |^T], 



^t(Q) = ^ 2 - \n t {Q) , n e (Q) = a 9 2 n £ a -a 



2 pq 

P 2~~ 



with 



m 



2 



g 2 NT 2 1 _ 1 1 

9 ' Ao= i^ ' A ° = (Q- P) 2 ' ?^^? T ? (9) 

For more details on the iT's (especially in our notation) see Appendix B of p3[ . If V — > oo , 
J2p turns into fd 3 p(2ii)~ 3 T J2n ■ We work with the Matsubara contour: Q = (iu n ,q), 
u n = 2-nnT . The gauge fields are Fourier transformed as 

A,{x)=Y.e- lPx A,{P) , A ll {P) = fe iP *A»{x) (10) 

with x = {—ir, r) and y = /(fdr J^ 3 ?"- To e.g. check this, the thermal Kronecker symbol 

P e i(Q-P)* = p VS nQ , np 5 q , p = [Q-P] (11) 



is very convenient. In (|]) the bullet on a. had been "forgotten", because the Greens 
function (|j) will turn out to be that of the trial theory as well. 

Were there not the paper ||, a suitable trial Lagrangian could come into mind while 
contemplating on (^). Use n free photon Lagrangians (numbered by a), supply them with 
variable mass terms such that @ is among their propagators, and identify the fields : A 
in @ = A in (0) : 

C = -\ F 2 + I A (MA) - (OA) 2 + icd 2 c (12) 
4 2 2a 



with = d^A v — d v A^. The trivial index a is suppressed here, and 

{MAY (x) = f Y^e- lQ{x - x,) [M t (Q)FT(Q) + M e (Q)W(Q)] A v {x') . (13) 

Jx ' Q 

The propagator of (]T2|) is (f|). But note that we are still free to choose e.g. constant 
masses 

M t = m\ , M e = m 2 ("m-case") , (14) 
or to cover the true leading-order propagators (|J) with 

M t = \ 2 t n t (Q) , M e = \j,n e (Q) ("A-case") . (15) 

Our trial Lagrangian (|1^) is non-interacting and quadratic in the fields A. The gauge- 
fixing term is necessary, because the mass terms are Abelian gauge invariant. To see this, 
insert the Fourier transform (PH|) into fllTf ) and notice that the gauge variation 8A V (Q) = 
~QvX{Q) drops out due to (|7|). Despite these neat properties of the mass term, the 
longitudinal one makes trouble. By — > id 11 in flllf ), and in the m-case for simplicity, 
we may rewrite A (MA) as 

A(MA) = A^x) \m\ 9^ v {%&) + m\ B^{id) } A v {x) , (16) 



In passing, (0) is the Abelian (trivial) case of the Lagrangian considered in [p4 ]. The 



matrix B (not R) has a denominator containing Q : V 2 = —q 2 Q 2 — > —AD. Hence, 
our trial Lagrangian has arbitrarily high powers in the time derivative. The definition of 
field momentum densities in higher derivative Lagrangians is a delicate matter, as is the 
construction of its Hamiltonian. Thus H, the trial object, makes the problem. Note that, 
if working with a constant-mass Stueckelberg term ^m 2 A fJj A fJi [25], this problem would 



not yet arise. We leave these difficulties right now, because there is a wonderful way out 
as detailed in the following subsection. 

II B. Feynman — Jensen 

To each Lagrangian, (|3]) and (0), there is a partition function which, using functional 
integrals, is expressed by the actions S, = — J 13 C, and S = — J 13 C : 

Z. = N (VA e~ s ' , Z = J- N tvA e- s . (17) 
Z' B J Zb J 

In (|TTD, and for the moment, let / T>A include the ghost field integrations. ( |T7| ) holds true 



in Euclidean space [^6], |21| : A is a purely imaginary field. The prefactors Z' B and Z B , 
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occur through the derivation of fllTD while integrating over 
8{dA — B) with normalized weight. Usually, they are hidden in the functional measure 
Af. But here, the two Af in (|T7|) are equal and independent of a. Nevertheless, they 
depend on (3 p6| . It might be emphasized, that the above functional language can still 
be applied to the Hamiltonian version V[H ] , since the first term of (jl|) is —Tin (Z), and 
the individual terms in (H, — H) (if known) can be related to Greens functions, which 
in turn derive from Z (with source terms included). 



Hamiltonians can be avoided at all, as we learn in § 3.4 of Feynman's text book |10 
(see also § 8.3,4 there and [[!]]]). Start from Z 9 , add the factors e s and e~ s under the 
integral, divide by (and multiply with) Z and define the average (. . .) as given in 
Then, 



Z. = !*<e-<*->) > |^-<— , (...)^^— • (18) 



The above inequality is, in the case at hand, the Jensen inequality, see e.g. p7| . Its 



simplest version states that (e ^ > e ^ . The proof rests on the convexity of e x [||, 12 
or, equivalently, on the nice figure 3.5 of [0. For convenience we take the logarithm of 

V [S] =F + T(S.-S)-T\n(^j =min. , (19) 

and call (|i~9|) the Feynman-Jensen variational principle. The last term, we come back to 
shortly, is obviously specific to gauge theory. 

In the non-Abelian case, there is a terrible pitfall hidden in (0). Admittedly, things 
were written down, to run into it with ease. If one still reads / T>A to include the ghost 
field integrations (wrong case), the ghost term would appear in the average T {S, — S) 
and loose the term linear in A, hence all A-dependence. But if the integration over ghosts 
is correctly recognized to be the Faddeev-Popov determinant, FP(A), it may be included 
into S, as S, nogllosts — ln(FP(A)). Now, in the avergage T {— ln(FP)), even powers of A 
survive. For the explicit formulation of this see Sec. V. We learn that ice, though being 
Hermitean, must not be viewed as a real number. Hence, in the non-Abelian case one 
needs to write the action S, in (0) and fllTf ) as to include FP(A), while there will be no 
ghost field integrations in / T>A. 

There is also the Peierls' version [28|, 0, ^, |5j of a thermal variational principle. It states 
that Tre _/3H > J2 n e~ <n ' H ' n> , rests on Jensen and may be used for another derivation 
H of (D). Things are closely related. But we have no rigorous answer to the question, 
whether the two versions ([I]) and fllTf ) are identical statements - just formulated in different 
language - or not (see also pt.3 in Sec. VII). The Feynman-Jensen variational principle 
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stays useful even at zero temperature. It has been applied at T = to A0 4 theory without 
Pm and with gauge sector |]18f . Formerly, these efforts had a Hamiltonian formulation, 
considered even thermal At finite temperature, but without the exotic last term, (|T9| ) 
has been recently used to obtain gap equations in lattice 4 -theory [13|. Before, it played 
a central role in a study of spin models and lattice gauge theory M . 



The role played by the unusual last term in 
gauge-fixing terms contained in S, and S : 



clears up by combining it with the 



V 



gauge 



l(^)/v>v-(nyi 



a 



- 1 - In 



ft. 



a 
a. 



V 



(20) 



For the logarithm in the first line see ( |A.14j ), the prime excludes n = q = 0. To understand 
the last term in the second line, remember that J2q = {T /V) J2 n S q - F° r the other terms 
insert (|10|) and use 

(21) 



{A lt (Q)A v (P))= [Q + P] G^{Q)= [Q + P] (RA t + BA^ + aDA c 



where the shorthand notation should be obvious from @. Now consider a of the trial 
Lagrangian to be one of the variational parameters. Clearly, with respect to a, (120) has 



an extremum at a = a., and V gange vanishes at this position. Moreover, it is a minimum, 
since the blank sum at the end in fl20 ) is positive (though quartic divergent). By far the 
best a is a.. We note three consequences of a = a,. First, the three terms selected in ( |20| ) 
may be simultaneously omitted in the sequel. Second, there is still dependence on the 
now common a, as it enters through (. . .) when traced back to the trial propagator (U). 
Hence, the above selection of a-dependent terms was incomplete. But the divergence of 
the last factor in ( pop helps maintaining the conclusion with rigour. Third, with respect to 
a, the variational principle is exhausted, so one should no more think about an "optimal" 
(common) a . 



II C. The "even version" 



So far, we were able to circumvent the Hamiltonian dilemma noted at the end of Sec. II 
A. But in the new version (|19D there is again a troubling element, as we become aware 
of next. Terms odd in the gauge field A (the AAA part of £. in particluar) drop out in 
V entirely, because they only enter (S, — S) and vanish there, since the average weight 
is the quadratic trial action. It is as if the 3-vertex was taken out from the outset. But 
a Yang-Mills theory with no 3-vertex can never be tested suitably by any trial-theory. 
For more details see the next Section. 
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For the resolution to this puzzle, it appears that the usual philosophy ("improve the 
trial theory") fails. Also, our trial theory fll2] ) is physically so reasonable: it "must" 
work. Our way out is to introduce one more version of the variational functional. On 
one hand this construction, which we call the "even version", is the decisive success in 
treating gauge fields variationally. On the other hand the idea is rather simple : in general, 
odd-in-A terms in the action can be avoided from the outset by playing around with the 
functional integrations over A as follows. 

Let us split the action into S, = £ + with £ keeping and O changing sign under A — > 
—A. The same decomposition can be done with the exponentiated action as e~~ £ e~° = 
e~ £ ch((D) — e~ £ sh(0). Since the second term drops out under the functional integrations 
over the gauge field A, we may write 

JvAe~ s ' = JvAe- £ + ^ ch ^ . (22) 

The new exponent, which we call —S„ , is an even functional of A. Since the above steps 
precede the use of Jensen's inequality, quite a new functional V arises : 

V [S] = F + T(S„ - S) = min. with S 1 .. = £ - In [ch (0) ] . (23) 



In (]23|) a = a, is understood, i.e. the logarithm of Z B 's is omitted together with the 



gauge fixing terms in £ and S. 

Once there are only even terms in the theory studied, the quadratic trial theory has 
a good chance to reproduce the leading-order perturbative results. We shall show in 
Sections IV and V that the "even version" works that way, indeed. There, the Faddeev- 
Popov determinant (depending on A in Sec. V, but not in Sec. IV) is part of S, and hence 
subject of the above "even"-ing procedure. 



III. TRIAL AND ERROR 



In this short Section we step back to the insufficient Feynman- Jensen formulation 
(Hp to see which way it goes wrong, to introduce some basic integrals and for a first 
run through the necessary algebra in the simplest case. For simplicity, let us even omit 
the Faddeev-Popov term (i.e. run into the pitfall noticed below (|19|)). It is not (solely) 
responsible for the defect, as we shall remark at the end of this Section. 



Using a = a, as reasoned below Q2TJ) , the functional reads V = F + T (S, — S). In 
the difference S, — S = — (C. — C) the terms odd in the gauge field A vanish under 
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the average (...). Others cancel. The only two surviving terms are 



T{S.-S) = V l -{A a {MA a )) + V 9 -f abc r s (AlA c v A^A vs ) = V M +V AAAA , (24) 

where y has reduced to f3V due to spacetime-independence of the averages. The first 
term, V m with M from (|13|) , is readily evaluated by using fllCf ), fl2l"D and the trace relation 
( [M t F\ + M e B] G )/ = 2M t A t + M e A e : 



V M = n \/T 4 (-L t - . (25) 

L t j are two sums out of the collection : 

Jt,t = -P 2 E \t > = E p2 (\e - A o) , *M = P 2 Y,P 2A ° A t,e ( 26 ) 
p p p 

with A tt £ = l/(-P 2 — Mt/(P)) and A = 1/P 2 . The prefactor n in (p5|) comes from the 
trivial sum over the colour index. 

The treatment of V aaaa starts with the Wick decomposition || of the average into 
three pairs with partners 

(Al{x)A b v {x)) = 5 ab Y,G^{P) = 5 ab T 2 \(u, v r + v, u s) . (27) 
p 6 

The first equality in (|27j ) derives with (|10D, (pip. The second one arises after integration 
over the directions of p. As the propagators A t ^ are rotationally invariant (even in the 
A-case (|l5l)), this angular integration amounts to the replacements R — > — |u, B — > 
— |w — ^p 2 A v and D — > ^(v — u) + ^p 2 A v with the Lorentz matrices u^ u , v^ u given by 
U^JJv — g^u and 4:U^U U — g^ u , respectively. For the sums r and s see ( p9|) below. Using the 
first equation (^7]) and with f abc f abc = nN one derives the first line of (28). Exploiting 
the w-f-version, one arrives at the second one : 



g 2 N 

V 'aaaa = nV —^~ 



E<V( p ) 



EG„(P)E^(Q) 



= nVT 4 ^ (r + s) (r - 2s) . (28) 
6 

The objects r, s in (P7|) , fl2"g|) are given by 

r = 2J t + J e + aJ , s = — - a J + aY , (29) 



where Jo and Yq are the sums of (|26[) taken at vanishing mass. 

The last term of V to be evaluated is the trial free energy F = — Tln(Z). First of all, 
since colours do not mix, Z is an n-fold product , 

F = 71 F co i our i css = -nT In (Colourless) = nV T A (-2I t - I e + I ) , (30) 
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and the colourless partition function is identical with that of scalar ED, see Sec. IV, if 
omitting the factor due to the scalars. In the formula ( |TT| ) for Z (read colourless and 
Euclidean), there are three unknown flying objects: A/", / T>A and / T>B. This is not a 
shame if Z is used exclusively as a generating functional. But here we need Z as a precise 
number. The trouble [^J with the normalization factor Af is proportional to the care of 
its treatment. We make efforts in Appendix A to write down at least (if not to derive) 
this factor Af '. Here we take from (|A.20|) that F indeed splits up into the terms in the 
right-hand side of (|3T)D. With ( |A.3j ), we obtain 
1 



2VT 3 



£ ln(-T 2 A,(P = 0,p))+E'ln(Po 2 A^)) , (31) 

P L n 

where the prime excludes n = 0. The index I may be replaced by t or by (then referring 
to zero mass). The expression (]31|) sticks with this awkward form as long as the A-case 



(15|) is included. But by differentiation with respect to A^ we may write 

- Xtdxjt = Li or h = h- f Xl d\ \ L £ (Xt = A) , (32) 

JO A 

the right half being equivalent to a coupling constant integration. In the m-case the 
above relation reads —md m I = L. 

The sums J, J to Y are divergent, and one has to keep track of variational-parameter 
dependences while renormalizing |J. To study this in simple terms (and for the rest of 
this Section) we turn to constant masses by (|14|) . In this case the frequency sum in ([31] ) 
can be done Using ( |A.7 ) and going to the infinite volume limit, one obtains 

e -y/x 2 +£ 



2tt 2 

Furthermore, Lp 



x 2 + e 2 e + In 1 



dxx 2 

ejJf . The sum J becomes 



ei = (3m e 



2vr 2 



dx 



x 



x 2 + ei 



+ 



1 J 



(33) 



(34) 



with clearly the |-term being UV-divergent as in (p3|). Even after subtracting zero-point 
energies by hand (which the functional integral does not know of), Ie — > h + -^z J^dxx 3 = 
J| ub , there remains a singular integral depending on the variational parameter eg. On the 
other hand, in a low-order perturbative treatment, such terms can be addressed as zero- 
temperature renormalization [|1~6| , p3f| and omitted entirely. As we like to reproduce these 
results only, the omission should be allowed here as well. Consider, for example, the 
combination — If nh — which occurs in F + V m, and supply p with an UV-cutoff A : 



rsub j 



-term 





,A/T 


47T 2 J 




/ dx 


4 








32tt 2 


\Te e 



1 x 2 + ei — x 



O e 



x 2 + ei 



(35) 
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Since we expect et ~ g , such terms are irrelevant in V up to g 3 . In the sequel we shall 
trust in the above arguments and omit the |-terms entirely. 



Deleting the divergent pieces this way (in e.g. (|33|) and ([34])), /, J, F become well 
defined integrals whose asymptotic series are known 



/ = 


IT 2 

90 


e 2 
-24 + 


J = 


1 

12 


e 

An 


Y = 


1 

8 " 


e 

-1T + - 

47T 



+ 



Me) 



ce 
64^ 



O e 



-dj 



(36) 
(37) 
(38) 



with e one of e t / = flint i, c — | + 2 ln(47r) — 27 and 7 the Euler constant. In the massless 



limit, the free energy 
radiation. 



is now recognized to be n times that of ordinary blackbody 



The contributions to V are now added up as F+ V m + V aaaa and filled with details : 

1 



V 



nVT 



-2I t -efJ t + 



—2L — e t Jt 



+ 



const + 



6 

nVT* ( 



a 



2J t + J e -Y e + -) (2J t + Je + 2Y k 



An 



£ -t-e t g 2 N 5 -±^ 
3 y 24 



1 + a 



+ ... 



(39) 
(40) 



There it is, the announced wrong result : V depends on a. Nevertheless, the structure is 
appealing : the parameters e t and £7 do not mix, the only extremum is a minimum, and 
its position has the right order g 2 N of magnitude. But, apart from this, the minimum 
positions e 2 t = g 2 N(5 + a)/2A and e\ = g 2 N(l + a)/8 give no sense : which a ? Including 
the FP-term, with the means worked out in Sec. V, does not help out of this dilemma, 
because it only leads to minor changes. To be specific, in ( |4~0|) 5 + a becomes 6 + a and 
1 + a turns into (2 + 3a)/3. 



IV. SCALAR ELECTRODYNAMICS 



For a first application of the "even" functional (|23|), we appreciate scalar ED as a 
suitable example. Remember that this system is an ideal toy model [I9~| to the gluon 
plasma, with view to the identical diagram structure, the need of resummation as well as 
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to its physical gross features. The Lagrangian, to be studied, is given by 

£. = (D^YD^ ~\f 2 -±- {dAf (41) 

with _D M = <9 M — igA^ and = d^A v — d v A v . By again identifying the fields (here : 
and A), the trial Lagrangian reads 

C = {dyyd,<P-- A F 2 -^-{dA) 2 + l -A{MA)-m 2 s <P*<P . (42) 

Its propagators are (£|) and 1/ [m 2 s — Q 2 ] for photons and scalars, respectively Here we 
concentrate on the spectum of real excitations. Hence, the mass matrix M is that of the 



the A-case (15). The variational parameters in the above trial theory are At, A^ and the 
scalar mass m s . The Lagrangian ([4~2"D turns into the effective Lagrangian (at order g 2 ) of 
hot scalar ED [TJj at the values At = A^ = 1 and m 2 = g 2 T 2 /A. So, within 0(g 2 ), the 



parameter space includes the exact answer (to be derived by variation). Note that both, 
original and trial theory, are invariant under regauging the photon field by 5 A = —dA. By 
definition, the decoupling ghost terms are kept apart from the above Lagrangians. But 
the Faddeev-Popov compensation must be taken into account in the partition function 
either by ghosts or as a determinant. 



IV A. The "even" functional of scalar ED 

Recalling Sec. II C, the partition function of scalar ED may be written as 

Z. = det\(3 2 d 2 ) Af I V {A^*4>} e~ s " . (43) 

The prime on the Faddeev-Popov determinant excludes the zero-eigenvalue (see also 
Appendix A). To specify £ and O in (p3|), in the case at hand, we read off from (|4l"|) that 

= -J P d with A = igAy*d^ - ig(d^*)Ay , (44) 

while — J £2 with £2 = g 2 A tJ- A fl <f)*<j) is part of 8 together with the quadratic terms in 
(f4l|). The index on £ refers to ^-powers. With a = a 9 , as required for ( p3|) to be valid, 
we may thus write the "even" functional as 

V = F + T^ l -A{MA) -m 2 s J%*(J)- J 1 " £ 2 -\n ch (j P a) \ , (45) 

where F = —Thx(Z), and the trial partition function Z is given by ( f4"3"D with all bullets 
stripped off there. 
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There is a high (but probably inevitable) price to be paid for the physical consistency 
reached with the above formulation : it obviously contains fields in arbitrary high powers 
(instead of only quartic). For the explicit evaluation of (f45|) one is, apparently forced 
again into a perturbative expansion, namely that of the logarithmic term ( In [ch(x) ] = 
x 2 /2— x 4 /12+x 6 /45— . . . ). But note, at least, that this expansion looks much simpler than 
diagrammatic thermodynamics : here the seagull vertex does not occur in higher powers. 
If, by any reason, terms of order g A may be neglected, then the functional simplifies to 

^trunc = F + T/f \a{MA) - m\ f 0*0 -f C *-Uf C i) ) • ( 46 ) 

In the following, while demonstrating the value of (fHf) , we shall in fact restrict to the 
truncated version ([46]). 



IV B. Evaluation of V trunc 

Let us group the above five contributions into "bare" and interaction terms: 

^ trunc = ^0 + Vint With V —F+ V M + V m 2 , V int = V A A4>4> + ^square • (47) 



For Vq we are well prepared from Sec. Ill : strip off the colour factor n from fl25|) and 
(p0|). Of course, with view to (|A.20|) , the free energy VT A {— 2I S ) of the scalars has to be 
added now. As the scalars have constant mass, I s is given by (^) with index i replaced 
by s. For V m 2 note, that the average (0*0) equals J2p S(P) with the scalar propagator 
given by S(P) = — 1/ (P 2 — m 2 s ) = — A s . Thus, in particular 

V m 2 = -Vm 2 s (0*0) = -VT% , L s = e 2 s J s , e s = (3m s , J S = -/5 2 ^A S , (48) 

and in total 

V = VT* (-2I t -L t -h- h e +I - 2I S - L s ^j . (49) 
Among the interaction terms, one is pretty simple : 

V AAH = -Vg 2 (0*0) = VT* g 2 (2J t + J e ) J s + VT* g 2 aJ J s = F™ . (50) 

Of course, Jo = 1/12 even in the A-case. To the right in d50|) , we have noted, that 
V AA<t><t> precisely equals the perturbative free energy contribution from the diagram co (one 
loop scalar, one photonic). But here the lines represent massive propagators, making oo 
depending on variational parameters. One may speculate, that the remaining term V squarc 
could correspond to the diagram (the inner line photonic). This is indeed the case, see 
(51) below. Classes of diagrams were whisked away in the treatment of Sec. III. 
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The first steps in treating V squave are straightforward : Fourier transform all fields, 
Wick decompose (the o— o-diagram drops out due to odd summand) and use fl2"T|). One 
obtains 



v square = -V \ g 2 E ( 2P - QY G AQ) (2P - QY a s a 

1 Q P 



(51) 

There is gauge-fixing dependence in V squSuTe through 



with A; = 1/ [ (Q - P) 2 - ml 

the propagator G from @. After some algebra one obtains that 



square 



n, (a=0) 
square 



-VT 4 g 2 aJ Q J s , 



(52) 



which cancels the a-term of Q5"0p. Thus, in the case of scalar ED and in its truncated 
"even" functional, there is no gauge-fixing dependence. This is true for all values of our 
variational parameters and for any mass matrix version. 

In ( 51 ) at a = terms with or Q v are projected out. The remaining sandwiches 
are PRP = - [p 2 - (pq) 2 /g 2 ] and PBP = P 2 - (PQf/Q 2 - PRP. Expecting the 
structure fl40|) , at least the terms linear in X t , X e , m s must be detailed. We therefore 
form differences as e.g. A t e = A t — A^, add ( |5T| ) to (pH|), write J t o = Jt — Jo etc. and 
split V int into a constant (v ), terms linear in such differences (vi) and the rest (v 2 ), 
which is certainly of higher order. Then, some terms of v\ (easily identified in (|56|) below) 
are regrouped into t>2, because they are of higher order by other reasons. Note that 
Ep (4P 2 Q 2 - 4(PQ) 2 + g 4 ) / [ Q 2 P 2 (Q - P) 2 } = 2 £p 1/P 2 . We obtain : 



Vint 

>'2 



Vg 2 (v + v 1 + v 2 ) 



with 
Q2p2 _ ( Q py 



Q 2 (Q-P) 2 J 

T 4 J (2J t0 + 3J s0 ) + 2 £ Ate{Q) ^ 



5T 4 
288 

1 



P 2 (Q - P) 2 



(pq) 



(53) 
(54) 

(55) 



2A t0 A sQ + 2 

Q P V 

(PQ) 2 



p" 



(pq) 



A t , (A S A; - A A ) + A, A s0 



P 1 



Q- 



A £ (A S A S 



A A ) + ^A, g 2 A A 



(56) 



In each term of (|56"D the argument of the first propagator is Q, and P that of the second. 
To understand why the last two terms of ( (5"6] ) are less than 0(g), the first one can be 
rewritten as —2m 2 J2q,p A^Ag A~ plus some products of differences. But for the last 
term in (|56D only a detailed analysis (of the type done in Appendix B) reveals its order 
T 4 g 2 ln(g) of magnitude. Such terms are known to occur in the perturbation expansion 
of the free energy |JI . 
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Up to order g 3 in V (or order g in Vx) only the line (|55|) needs further study. Note 
at first that, formally, the expression (§) (at N — 1 here) appears in this line. So, the 
machinery "knows" of the leading-order longitudinal polarization function. There are 
two ways to evaluate Y^o A^ t ne. (for later use, we detail both). First, one may cancel 
i7-functions with those in the trial propagators (using IIi = 3m 2 — 2i7 t ), write v\ in 
terms of basic integrals as 

£ A et n e = — T 4 ( ^L t + ±L e - 9 jJ t ) (57) 

and proceed with expanding the latter (see below). Note that these cancellations are 
possible in the A-case only. The alternative second way is by far the easier and more 



enlightening one. As is basic to the dimensional reduction method p2l Q and to 



various related thermodynamic calculations (e.g. [JTJ), a frequency sum may be 



occasionally reduced to its Qo = term. This step, if valid [|T^, |56fl , rests on the structure 
of a massive propagator A(Q) and usually prepares its soft part while contributions from 
non-zero hard frequencies are of higher order. Of course, there must be no hard part in 
a sum under such study. In fact, ([57]) is an ideal example for the above. Moreover, at 
Qo = the polarization functions il^ and H t reduce to constants, namely 3m 2 and 0, 
respectively : 

Dependence on X t has dropped out, with the reason readily detected in the vanishing 
factor II t of A 2 , i.e. in the absence of a (squared) magnetic mass at order g 2 . In the same 
manner the soft parts of J-integrals are obtained : 

J e = Jo-^^ + ... , J t = Jo + + • • • (59) 

with again no dependence on X t by the same reason. 

To complete the evaluation, note that the scalar contributions J s , L s , I s are m-case 
objects, hence their expansions are given by (|36|) and (|37D with e = e s = (3m s . For the 
remaining A-case integrals in V , (f4D]), apparently, the sums L t and are still to be 
studied separately. For this somewhat delicate task see Appendix B. As a result, both L 
start with a g 2 A 2 -term whose prefactor k diverges logarithmically. From ( B.5|) to (|B.7|) : 



1 3 47^^/ V24 J 3 V ; 

Now note that the singular piece k drops out in V , because there L appears in the 
combination 2I+L = 2/ +|l — 2 J^dX j\ L, see (|32D, and the curly bracket is a projector : 
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{ } A 2 = 0. One might ask for the fate of the singular terms in fl57p . They drop out there 
by cancellation, and Q58| ) derives again. 



IV C Minimizing V tTunc at order g 3 



In the preceding subsection, the expansions were driven just as far as to allow for 
writing down the functional V up to third order in the coupling g. Of course, as in Sec. Ill, 
we anticipate that the solutions to m s and will be O(gT) and 0(1) in magnitude, 
respectively. By combining the details of the preceding subsection one obtains 



V 



to g A 



I 45 288 24tiV3 



\ 3 

— - \ e 
3 



+ 



1 



12tt 



3o 2 

(f3m s f - -^(3m s 



(61) 

It still has the structure of fl4*0|). The variational parameters do not couple, which is specific 
to the order considered. The absence of any dependence on A t was already understood, 
although merely technically (see also Sec. VI). The above V, when plotted over the A^- 
A t -plane, has the form of a long gutter. The resolution of this defect is deferred to Sec. IV 
D. 

Minimizing ( pTj ) with respect to m s and Xe gives the values 

1 



gT 



A5 



(62) 



as expected. We immediately look for the value of the above V taken at these parameters, 
which is the height of the bottom of the gutter : 



V 



VT 4 



71 



5g 2 

-2 — + — 

45 288 



1L 
12n 



3V3* 



(63) 



with the last term in the square bracket being due to the scalars. The minimum perfectly 
agrees with the perturbative free energy up to g 3 . The g 3 -term, the correlation energy, 
was given by Kalashnikov and Klimov |]37| (eq. (19) there, taken at A = /x = and 
e = g). In summary, for scalar electrodynamics and up to the third g power, the "even" 
variational functional has all required properties, namely gauge fixing independence, the 
right minimal value and (apart from degeneracy) the right minimum position. 



IV D. Solution to the gutter problem 



The missing dependence on A t in (BT) is, as already noticed, an artifact of the restric- 
tion to order g 3 of the functional. The problem merely is how to go one order higher within 
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the expansions so far developed. First of all, we notice that g 4 -terms are allowed within 
the truncated functional, although the neglected next term of In [ch(a;) ] does contribute 
at order g 4 too. However, the latter is a constant at this order; variational parameters 
appear at g 5 . 

Let us try to avoid expansions, and let the collection of all terms containing At be 
denoted by V^l nc . Up to an additive constant, it may be written as 

V£L = y^U t + Vg 2 vf with U t = -L t + 2 f^dX ~L t (A t = A) + ^L t . (64) 

Jo X Af 

Here, ?4 is made up of the first two terms in (|5lJ), but leave aside for a moment. 
Then, the minimum condition may be given the form of a product 

= d Xt Ut = - dx t Lt ) ■ ( 1 " T2 ) ( f > r ) ) 




vAt 

with the first factor "unknown", but the second reaching zero at X 2 = 1 as desired. To 
be sure that this zero corresponds to a minimum, the first factor must be shown to be 
positive. We shall do so at the end of Appendix B. There, the first factor is also seen 
to be of order g 4 and to vary as A^ for small g, see ( |B.11|) . Hence, U t has a Higgs-type 
shape Ut ~ const — g 4 X\ + \g^X\ + . . . with a maximum at the origin. The curvature of 
the gutter sets in one order higher, indeed. A plot of V now merely looks like a long 
bath-tub. 

The above construction only works if the correction v% remains below the order 0(g 2 ). 
Its first term is the first in (BB) and is of order T 4 g(J t — Jo) in magnitude. With view to fl55| ) 
it is indeed below g 2 . For the second contribution (the second in ( |56"D but with A t0 in place 
of An) we need a bit of calculation. Both sums may be considered "soft" , i.e. n(x) — ► T/x 
is allowed, thereby preparing the contribution of interest. All propagators are represented 
spectrally. For the two frequency sums, eq. (6.6) of p3) is used repeatedly. The result is 
a 3-momentum double integral over (among other factors) Jdx^ (^pi X \x,q) — p^°\x,q)^j. 
But, due to the sum rule ( |C6| ), this factor vanishes, q.e.d. 



V. YANG-MILLS FIELDS ( the gluon plasma ) 



For treating the non-Abelian theory (|3]) in its "even version", we use Sec. IV as a 
guideline. Hence, first of all, we strip off the ghost terms from £ and introduce the 
index " no" for such reduced Lagrangians : 

S. no = - fc™ = - f{Co + Ci + C 2 ) with (66) 
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d = -g(d fl Al)f abc A^ h A vc , C 2 = ~g 2 f obc r s A b lx A c v A lxr A vs . (67) 

Here, £ is the quadratic part of hence including the gauge fixing : C = — (F a ) 2 /4 — 
(dA a ) 2 /(2a,) . The Faddeev-Popov determinant now depends on the gauge field and is 
thus subject to functional integrations. But for convenience we may split off its bare 
factor. The partition function, still waiting for its even-odd decomposition, so far reads 

Z. = A-det' ((3 2 d 2 5 ab ) M J VAle- s -°- SFP . (68) 

The two factors in (0), which obviously stand for the FP-determinant det' (fl 2 dD), derive 
through 

det'(p 2 dD) = det' (p 2 d 2 5 ab )det' ( \ d 2 5 ab - d» g f abc A c " ' 



d 2 5 ab , 

= det' (f3 2 d 2 5 ab ) det' (1 + W) = det' (f3 2 d 2 5 ab ) e~ SFP , (69) 



where by W the part odd in the gauge field is prepared : 

1 



W = -gf abc d»A c - 2 . (70) 
The first acts on A c ^ and all functions that follow. We read off from (^) that 
S FP = -In [det'(l + W)] = -Tr'ln(l + >V) 

= -5 Tr '" 1 ( 1 - w2 )4 Tr ' 1 "(r ± w) • (71) 

In the second line, clearly, the even-odd decomposition is achieved. But the second 
equality in ([H]) (first line) is delicate, because all eigenvalues of 1 + W have to be positive, 
but are not. While this point needs care in exactly solvable models [j38H , here we may be 
content with a crude argument. For the intended comparison with perturbation theory, 
the above logarithms are expanded anyways. Hence (|71|) is merely a formal compact 
notation for series to be generated [|39| . 



We are ready to form the non-Abelian "even" action S„ through + Sfp — > S„ = 
£ - In [ ch (O) ] with £ , O given by 

S = -|Vo + A)-iTr'ln(l-W 2 ) , (72) 

rP 1 /1 + W\ 
O = -/ A-^'ln^) . (73) 

The trial theory has remained unchanged. It is that of Sec. III. The trial partition function 
is given by (H) without the bullets, at S FP = and with S no = - f C no . The free energy 
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F is (0). Thus, the "even" functional fl23|) of the gluon system (taken at a = a,) reads 
V = F + TU -A a (MA a ) — J C 2 - ~Tr' In (l - W 2 



In 



ch 



1 



r l + W 


>] 


Li-w. 





(74) 



The first two terms form the bare part an d are familar from Sections III, IV: 



V = nVT A (-2I t -L t -h-^L e + I ^ • (75) 

As in Sec. IV, we expand the logarithms up to W 2 to reach a reasonable simple 
"truncated version". Since Tr'W = 0, no such term arises from the last logarithm. 
Thus, 

^trunc = ^0 + Vint , Vint = V AAAA + V FP + v 

square 

(76) 

with 

V FP = ^(Tr'W 2 ) , V^ = ~((fc i y j S j • (77) 

The contribution V aaaa is given by (^). It agrees with the perturbative free energy 
contribution from the tadpole diagram (both lines gluons) : V aaaa = F™ . Compared to 
Sec. Ill, there are two additional terms in ( f76| ) : the last two. By analogy with Sec. IV 
we expect that they equal the two other diagrams at second order, which were missing in 
Sec. III. Indeed, taking the trace of W 2 with states (f3V)~ 1 ^ 2 e~ tPx , using f abc f abc = N n 
and through Wick decomposition, we obtain 

v FP = ^yg^MQ) ^:^ = K , (78) 



v„ e = nv y —Y J Y.[(Q + p ) XG UQ-P)G up (Q) 

1 Q P 

-G Xt {Q)G Xt {Q-P)Q p ] G m {P){2Q-Py = Ff , (79) 



where, in (|78|) the symbol " V (with two out of many dots) stands for the ghost loop with 
an inner gluon line. 

Quite different from scalar ED, the gauge-fixing dependence does not cancel in a 
manner independent of variational parameters. Splitting the Greens function as G = 
X + ctDA , we see that a occurs up to the third power. The term a 3 is contained in 
^square only, and its prefactor vanishes. Collecting a 2 - and a-terms one obtains 

V (a2) = -nV 9 ^a 2 EE p4 / P) 4 P VWr - (the same at zero mass) , (80) 
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+ Q 2 (P 2 - Q 2 ) X^(Q)X^ (P) ] ~ (the same at zero mass) . (81) 



The fact that ( jSCf) and (pl|) vanish at zero mass reflects gauge invariance of thermodynamic 
perturbation theory at order g 2 . For the next step, namely analysing Vtmnc at order g 3 , 
we need more: fl8~0|) and (81) must remain below g 3 . This is the case, as one may check 



e.g. by power counting. Remember that perturbatively a g 3 only arises by dressing the 
g 2 -diagrams, whereby gauge invariance persists. 

The strategy of further evaluation is now that of Sec. IV, as detailed above (0). Since 
they are of higher order, we temporarily omit the two a-dependent terms (jSOf) and (|8lD . 
In Vi nt this amounts to the replacement G —>■ x = + BA^. Then the terms (t>i) linear 
in x ~ Xo are isolated, and terms of higher order - others than in Sec. IV - move to t>2. 
But evaluation of v i runs through the steps in Sec. IV B and, surprisingly, ends up with 
the same result as in Sec. IV, namely (0) at m s = 0. Just to show prefactors : 

rp4 rp4 rp4 rp4 

Vint = nVg 2 N (v + v x + v 2 ) , v = — , v x = + 2g 2 NX 'j Le ~ ~Q J ° ' (82) 

The complete functional up to order g 3 (add ( ^2|) to (|75|) ), does not depend on X t (gutter 
form) and reads 

y to 9 3 = const + nVT^Ue (83) 

with the function Ui defined as U t in (0) by changing the index. Minimization gives 
Xe = 1, as desired. For the height of the minimum to order g 3 we obtain 



n i min T rrri4 

U = nv 1 



45 + 144 ~ 127 I ) 



This is equation (8.47) in |^lj| . At N = 1, the correlation energy (g 3 -term) agrees with 
the photonic one in scalar ED, see (|63|). 



As in the Abelian case (Sec. IV D) the functional is expected to become convex with 
respect to A t by including g 4 -terms. However, at this point we run into non-Abelian 
difficulties. There are four terms to be included. The first one is U t (replace Ut in (|83| ) 
by Ui + 2U t ), which has a minimum at X t = 1. The second term arizes from t> 2 in (H), 
a rather lengthy expression (seven lines say) and so far not evaluated. The third and 



fourth terms are the a-dependent pieces (|30D, (|81]) and cause the trouble. They should 
be (but are not) either constant, or minimal at X t = 1, too, or of lower order in magnitude. 
Consider e.g. the A t0 -part of the a 2 -term (BO). If evaluated "soft" it vanishes (in the 



manner noted at the end of Sec. IV). At first glance, as no UV-cutoff is needed, one might 
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/ 2 \ 

conclude that V\ a — at all. However, it appears that there is still a hard contribution, 
which in turn needs no IR-cutoff. Because this is perhaps somewhat unusual, let us state 
the result : 



V t ^ = -nVT^ 9 -^±Z with Q= r dx ' I dt ^ 

*> O A A O — A /„ 1 /„ i. 



24 48vr 4 Jo e x - 1 Jo e \ xt - 1 

\wh + (^ + (TTTf la{t + 2)+ (^W la]t - 21 ) ' (85) 

The derivation (a mess) used ( U.2 ). To check the above statement of vanishing soft part, 



one may write 2/(xt) for the second Bose function. Then the integral over t gives zero, 
as required. But as it stands, S is some non-zero mathematical constant ~ —1.04). 

The above remaining a-dependence, which prevents us from solving the gutter prob- 
lem in the non-Abelian case, is the "minor detail" noted in pt. 8 of the Introduction. 
There must be a resolution to this puzzle within the truncated version ([F(]), because the 
terms beyond, depending on A, are of order g 5 . As the vicious term (^) contains two 
Bose functions, the way out has probably nothing to do with renormalizations. The only 
possibility we are able to invent is the fact that at higher orders there is also a C-term 
(see ((D) in the propagator, which is missing in (f|) and is specific to non-Abelian the- 
ory. Furthermore, this term has a factor a in front of it, see e.g. § 3 of [|0|. Let such 
speculations be beyond the scope of the present paper. 



VI. STATIC PROPERTIES 



So far, while testing the "even version" in the A-case, we were thinking in terms of real 
excitations in the plasma (scalar and gluon), whose spectra are hidden in the polarization 
functions. Here we recall the other well-tractable case within the infinity of Abelian gauge 
invariant mass terms. Before all, turning to the m-case comes with a change in philosophy. 
We now ask for the best constant-mass terms (longidutinal and transverse) in the trial 
Lagrangian. To leading order (otherwise see e.g. [0), static propagators have the form 
(— q 2 — ^screen) 1 - But the trial propagators read (Ql — q 2 — m 2 ^ 1 . Nevertheless, it may 
well happen (remember the n Qo = 0-method" of Sec. IV B) that they loose memory to 
their dynamical element Qq automatically. 

For Yang-Mills fields, the analysis runs through the steps of Sec. V up to (|82|) . No 
gauge-fixing dependence occurs up to the order g 3 to be considered here. The bare part 
is given by (|75D, now with the m-case integrals (|36|) , (|37D to be inserted. The crucial 
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line where the m-case starts to make differences reads 



v x = 2T 4 J ( J t - J„) - 2 £ Aft(Q) n e {Q) . (86) 

Q 

Within the present accuracy, the above sum may be reduced to its Qo = 0-term. But note 
the difference to the A-case. Once the transverse propagator is supplied with a non-zero 
magnetic mass by hand, this variational parameter survives in the result : 

= — ^— ((3m e - (3m t ) . (87) 
Q l2lx 

The same happens in the J-sum, see (p7|). But the combination of these details in (|86| ) 
yields v\ = — T 3 m^/(247r). The linear (not the cubic, see below) dependence on m t has 
gone, this time by cancellation - a wanted detail, as we see next. Including the bare part 
Vq the functional reads 

V t0 m 9 -Le = nVT 4 (-^ + 9 -^L + -L\(p m ) 3 -g 2 NPm e } + -^-(/3m t ) 3 ) . (88) 
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The longitudinal part clearly becomes minimal at mi = gy/WT/y/W , which is the well 
known Debye screening mass at leading order. There is a transverse part in fl8B|), hence 
no gutter problem. As m t is restricted to the positive half-axis, the minimum is reached 
at m t = 0, which is the magnetic mass at the order studied, indeed. 

In spite of the above correct answers on static properties, there remain delicate ques- 
tions. Remember that the (squared) Debye mass 3m 2 already entered the dynamical 
calculation at ( j58|) . It appears that, within the order g 3 , the variational functional can 
not really discriminate between statics and dynamics. In fact, the minimum value of the 
functional (|88|) agrees with fl8~4|), i.e. with the exact one to order g 3 . Thus, two equally 
low minima are found over the space of mass terms. However they are joined, namely 
through a subspace of all functions Hi that have the value 3m 2 at zero-frequency, and 
II t vanishing there. Nevertheless, in the A-case the appearence of constant masses is a 
technical byproduct, while in the present static case it answers the posed question. Let us 
add conjectures on the behaviour in higher orders. The safe ground is on the dynamical 
side. Supplying the variational functional with anything good then it might answer with 
self-energies comparable good. For static properties, on the other hand, one needs more, 
namely some philosophy of why the trial propagators get rid of its dynamical part Qq 
by only forcing the mass to be constant. Remember also that, starting from the real- 
excitation spectrum in the tu-g-plane, the static limit (u = 0) is only reached through 
a range with imaginary wavevector [SB] on mass-shell lines. Perhaps, the variational 



procedure prepares at least the first non-vanishing term of each screening mass. 
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At the supersoft scale, the magnetic mass (see |33| , fill for more recent work) most 
probably comes with some numerical factor times g 2 T 
the last term in 



Then, as a rough speculation, 



V = nVT 4 ^const + (/Jm^ 3 — const g 4 (3m t ^j 



(89) 



would be in search. Note that such a term, if any and if no others, would arise in one 
step over the present truncation of the functional. For possible danger with this step see 
the last point in the following list of open questions. 

For completeness, we add the m-case result for scalar electrodynamics. It simply 
agrees with (|88[) at N — n — 1, except for the constant terms and an additional term due 
to the scalars, which may be both read off from (j6T|). Let us end up with the question 
which way the magnetic sectors of Abelian ji3] and non-Abelian theories might become 
differenct in a variational treatment. 



VII OPEN QUESTIONS 



In the preceding Sections, the application of the variational calculus to pure gauge 
theory was far from being a straightforward procedure. Several problems were eluded and 
questions not answered, because we could not. Let us recall these questions and just list 
them here. 

1. The Hamiltonian formulation to both, the Gibbs-Bogoljubov or Feynman- Jensen 
varational principle (see text below fli7|)), was given up in Sec. II because we were 
unable to construct the Hamiltonian H of the trial theory. This construction is a 
challenging task. See the text below (|16|) . 

2. Knowing the Hamiltonians of both, trial and studied theory, one could construct 
the common physical Hilbert space. By forming the BRST-charge and projecting 
out physical states from the outset, this would be the natural approach to the 
Gibbs-Bogoljubov version || ^|, 

3. The functional V in both versions, Gibbs-Bogoljubov and Feynman- Jensen, has the 
total minimum value in common (namely the exact free energy). However, the trial 
spaces are different. Hence, a given trial theory which does not cover this minimum 
could lead to quite different approximations. Since presumedly, this is not true, a 
proof of the full equivalence of the two principles is desirable. Note that such a 
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proof would circumvent our Hamiltonian problem of the above point 1. Moreover, 
the interpretation of the trial space as one of non-equilibrium statistical operators 
would be preserved. 

4. We have not made an effort to introduce, by Legendre transformation, the 1PI- 



generating functional T, although there is a variational principle even to T [f44l E9 



Renormalization ||, not yet needed in this paper, is probably inevitable already 
when the method should reproduce the next-to leading order perturbative results, 
such as e.g. the lowering "by glue" of the longitudinal plasma frequency (for scalar 



ED this is the term — 0.37e in eq. (5.5) of fT9fl ) 



6. From subsections II B to II C we turned to the "even version" immediately. But 
perhaps there is something in between that we have not found, namely a feasible 
modified trial theory not running into the pitfall of Sec. III. 

7. Only a very poor subspace of polarization functions was considered by simply vary- 
ing prefactors X^t in front of the true functions n t , already known perturbatively. 
A honest "even version" -variational treatment might instead vary unknown func- 
tions II t i(Q). To make sense, this generalization probably needs g 4 -terms in the 
functional V . 

8. For the "minor detail" of reminescent a-dependence when solving the gutter prob- 
lem in the non-Abelian case see the comments at the end of Sec. V. 

9. The most terrifying step in Sections IV, V was the expansion of the In [ch( ) ] term 
in the variational functional. So, the question is whether this expansion can be 
avoided some way. 

10. With regard to the observed gauge-fixing independence, it could turn out that a 
later truncation of the series makes less sense than reading In [ch(x)] ~ \x l as some 
good approximation. 



VIII. CONCLUSIONS 



A Feynman- Jensen type thermal variational principle is constructed such that an 
Abelian free trial theory works well in both cases, scalar electrodynamics and pure Yang- 
Mills theory. To this end their actions are to be rewritten such that only even powers in 
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the gauge field appear. This way, the perturbatively known leading-order self-energies 
of photons, scalars and gluons, respectively, are reproduced (apart from a minor open 
question to the non-Abelian case) by variation of their prefactors. The subspace of 
constant masses covers the inverse Debye screening length. There is a large asymmetry 
of the functional with respect to the (photonic/gluonic) transverse sector, as it does not 
(yet) depend on the corresponding parameter at order g 3 . 

The delicate problem of handling two different covariant gauge-fixing parameters (one 
of the original and one of the trial theory) has a simple resolution : they become equal by 
minimization. Hence, the observed gauge-independence refers to the remaining gauge- 
fixing parameter common to both theories. 

The new variational functional contains a term In [ ch (AAA) ] and hence involves 
arbitrarily high even powers of the gauge fields A. In the non-Abelian case (and within 
covariant gauges) such powers occur already in the unmodified Feynman- Jensen principle 
due to the Faddeev-Popov determinant depending on A. Unfortunately, for evaluation 
and minimization we had to expand the ln-ch-f unction. But a true nonperturbative 
scheme should never refer to g-powers at all. So, the present success is still below the 
potential nonperturbative possibilities of the variational approach. 
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Appendix A 



Here the functional integral measure of the trial partition function Z of scalar ED is 
made explicit. Z is without the dots there. The normalization factor M is fixed by 
requiring that, in the massless limit, the partition function Z must turn into two times 
that of blackbody radiation, one of the photons and one of the scalars. On the more 
ambitious task of a true derivation see the comments at the end of this Appendix. 

We start by splitting Z into four factors, Z = Z a Z^et Z& Z s with a piece of M contained 
in each. But notice the redundance of such a factor in front of an unspecified JT> . . ., 
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hence e.g. Z a = 1/Zb suffices. The simplest part is Z s = "J T> {0*0}" e Ss with 



(m 2 0*0 + 0*^) = £ (m? - P 2 ) 4>(P)*4>(P) = £ „ • 

P p , 71 P " 



(A.I) 



At each of the countable infinite discrete positions p , n there are, as is complex, two 
independent integrations. ( |A.1|) refers to our convention 4>(x) = e~ lPx (f)(P) but we may 
turn to that of Kapusta [^T[ by (f>(P) = \/(3V'(a + ib)/^/ 7 ? (with indices p, n on a, b 
suppressed). We now guess the functional integral measure and justify by evaluation: 

Zs = MoU^p Jdadbe-^- p2 )( a2+b2 ) = ^l[(-T 2 A s ) , (A.2) 



p,n 



where T7 = Tip n an d 



Mo = n n' 2 ™ 



(A.3) 



p n 



with the prime excluding n = 0. Remember that A7 1 = P 2 — m 2 s = — (27mT) 2 — p 2 — m 2 < 
0. Of course, each factor in Ao has to be attached to the corresponding one in the n— 
product in (|A.2|) , and the product over n has to be performed first (other constructions 
may be possible). 

The infinite product ( |A.2j ) can be performed. To this end we collect four (known) 
formulas of general use. By contour integration : 



' E p2 ,.2 



n 1 



1 

X 



+ n(x) 



(A.4) 



with n(x) = 1/ UP* - 1J the Bose function, (gjg ) is eq. (2.38) of 0. Multiply (TAJ) 
with 2x and integrate over x from c to y : 



V _ p2 

y r o 



c 2 -P 2 



y-c + 2Tln 



1 - e - ^' 



1 - e-^ / 

Multiply ( |A.5| ) with /3, set [3y = uj and perform the limit /3c — >■ : 



(2nn) 2 J 

Exponentiating (|A.6|) and extending to all n, one arrives at the fourth formula 

• (2vm) 2 



(A.5) 



(A.6) 



n 



u 2 e- w 



u 2 + (2vm) 2 (1 - e - w ) 



2 > 



which is eq. (89.5.16) in @ and (2.269) in p[. Check 
( |A.2|) we obtain 

z s =um -t 2 a s =n 



(A.7) 

at u; — > 0. Using (|A.7| ) for 



i.e. 



(Ai 
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F s = -Tln(Z s )=2j2 



m 2 + p 2 + T In ( 1 — e 



(A.9) 



which is, at zero mass, the desired result of twice a half blackbody radiation. The guessing 
was good. Aside, one could include the zero-point energies by the redefinition Mo — > 
n p e^ 2 nn'2vrn. 



We turn to the factor Z dct with again a trial-and-error prefactor : 
Z det = M - 2 det> (fd 2 ) = [^n't-TV" 1 



L P 



nn 

p n 



/ {(3pf + (27m) 



(2vrn) ; 



(A.10) 
(All) 



where in the blank n' an d on the determinant the prime excludes only the one position 
n = p = 0. As the determinant is the product of the eigenvalues —(3 2 P 2 , naively, P = 
must be excluded to make sense. However, if this is required to result from a derivation, 
one might go back to the unity-insertion in the Faddeev-Popov procedure : 



1 = A ■ Mq 



tf T 2 da 



[IT* 



-P'a 



(A.12) 



Originally the ^-argument was d 2 A (with SA^ = —d^A the gauge variation). Since there 
is no constant term in A, there is no ao,o~integration in (|A.12|) and no zero P in (|A.10|) , 
q.e.d. (TA~T^ ) directly leads to A = Z det . Using ( pO) ) for (TA~Tl| ) we have 



det 



-T In (Z det ) = -2 Y! \ 5 P + T In (l - e 



-Pp 



(A.13) 



Clearly, with the above measure, the determinant-term subtracts twice a half blackbody 
radiation. In passing, the prime in ( |A.lij| ), while still being necessary there, becomes 
irrelevant in the continuum limit. 

With an argument quite similar to that below ( |A.12j ), there is also a prime in the 
measure of / T>B. This integration runs over a 5 (dA — B). But dA cannot be constant, 
since otherwise A would be linear in spacetime and lie outside our space of Fourier- 
transformable fields. So, P = may be excluded : 



VB " e"2 



±Y! P B{PTB(P) 



right 

n 



2tt 



da db e 2 ° 



IT 



(A. 14) 



As B(x) is a real field and B(-P)* = B(P), B{P) = v^F(a + ify/V?, the two inte- 
grations are placed on half of the P-space, the right say (let right and left exclude the 



origin). The prefactor was chosen here to reach the simple result Z a = 1/Zb = H'y^/ a 
It must wait to make sense in combination with Za- 
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The photonic part of the trial action includes the mass terms Mt,e '■ 

Sa = W{P 2 -M t )A-^A v + W{P 2 -M,)A-^A u 

£ p P 

+ ^E p2 ^ DM ^ = S a + S a (A. 15) 

with A~ = A^—P) and the D-term being part of S A . The corresponding further splitting 
Z A = Z A Z A is allowed because the transverse components (those in e 12 -direction, e 1 )2 _L 
p, ei _L e 2 ) in the expansion 

A"{P) = u x E^ +u 2 E% + vT» + iwU^ with T* 1 = ^0, ^ , E^ 2 = (0, e 1)2 ) , (A.16) 

drop out in S l A and are the only parts surviving under the R-operation : R^Ei^ v = E^ 2 . 
As the first three terms of ( |A.16| ) as well as wV^ are Fourier transforms of real fields, 
half P-spaces are related by Uj(-P)Ef(-P) = Uj {P)*E?(P), v(-P) = -v(P)* and 
w(—P) = w(P)*. Hence, the integrations in Z A , to start with, are of the real-field type 
( A.14| ), except that there are now two integrations at the origin n = p = and four at 
each place in the right half. Two of the latter may be attached with the left half. Then, 
choosing the same functional integral measure as for Z s , we arrive at precisely ( |A.2|) with 
the role of m 2 taken over by M t (P) : 

Z A = ^Y[(-T 2 A t ) . (A. 17) 

The longitudinal part of the action is first rewritten as 

* Ifi'-^W) - th ? = p T ! f . (A.i8) 

2pl/^\ L J a J rj = —iPqW + pv 

Next we observe that £(— P) = £(P)*, rj(—P) = —rj(P)*, and mark the origin and the 
right half P-space to count independent integrations (two over £ at the origin and four 
in the right). Finally, by changing the variables from v, w to £, rj (with unit Jacobian 
determinants), and with the now familiar functional integral measure, one arrives at 



4 = ^oII \I-T 2 ^ Ao nV- T2 « A o' • (A.19) 

Note that most of the above "trivialities" were due to carefully counting all positions in 
P-space, i.e. to place the primes right. 

We are ready to constitute the scalar ED partition function from the above several 
factors : 

Z = IW A/q IT (-T 2 Ao) 



K 2 I[(-T 2 ^) V[n«] A? I]' (-T 2 A ) A" 2 n(-T 2 A S ) . (A. 20 ) 
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Obviously, the gauge-fixing parameter a cancels. Now, counting halves of blackbody 
radiation amounts to —2 + 2 + 1 + 1 + 2 = 4 as required. 

A true derivation of the above must not anticipate the known zero-mass results. With 



|26[ as a guideline, such derivation should be possible even inside covariant gauges, i.e. 
without a recourse to physical gauges. There is one problem in taking the right starting 
point (maybe with a factor Mq in front of the classical partition function for the four [of 
six] degrees of freedom to be quantized), and in the volume factor (to be split off) the 
other. 



Appendix B 



Here the two sums L t and L t are evaluated, in the A-case and with regard to contribu- 
tions not accessible by a naive Qo = 0-method. The details are required for subsections IV 
B and IV D. We start from the definition (EBI) and work with the spectral representation 



1 



P 2 - \ 2 n t , e (P) 



P 2 - X 2 



(b.i; 



of trial propagators. The above spectral densities are related to ordinary ones, denoted 
by pu(x,p;m 2 ), by 

(B.2) 



P$( x iP) = Pt,£ (x,p;X 2 m 

Hence all sum rules ( |U.6| ), ( |C7| ) remain valid for p^ if m 2 is replaced by A 2 m 2 to the 
right. Using ( |B.1|) and the sum rule 1 = Jdxx p^ x \x,p) the L-sums read 



Lt,i = -/3 4 E ( p2 ^ - l) = -/? 4 £ Jdxxpl x J(x,p) 



' Pj-P 2 
P 2 -* 2 



- 1 



(B.3) 



Next, with ( |B.2| ), defining p = (x 2 — p 2 ) p (cf. ( |U.1|) ) and using (|A.4| ), we may write 



P 4 



2vr 2 Jo 



dpp l dx p te (x,p; A m ) 



+ n(x) 



(B.4) 



Note that both, p and the square bracket, are odd functions of x. Let us split p into its 
leading part as given by (|C.2|) , (|C.3|) , and the rest ~p — p lcad ; which we call the soft part. 
Correspondingly, L is written as L lead + L soft . Introducing an UV-cutoff A, the leading 
parts may be written as 



^lead 
^lead 



A? 



2 n 2 M 



P 2 9 



2X 2 12vr 2 Jo 



1 

- + n(x) 



1 - I(x) 
I(x) 



Kx) 



dp - . 

P 



(B.5) 
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(P-5|) is the right place deleting the "|-term" as discussed below (|34D in the main text. 
But even under the control of the Bose function there remains a logarithmic divergent 
factor, namely 

P 2 



dxxn(x) 



Using 
in (|6 



we have L' ead 



47T 2 

- k) \g 2 X 2 and L^ ad 



f dp- . 

Jo p 



(B.6) 



Kg Af, which are the A -terms 



We turn to the soft parts of L t} £, whose series might start with g 2 X 3 . To prepare this 
A 3 -term, one may simply write T/x in place of the square bracket in ( p.4j ) (and, of course, 
the mentioned difference in place of p). Using the sum rule (C.5), one obtains 

-I i / ^ \ 3 

-£° /"CO fir "1 

/ dpp 2 dx- p e (x,p) -p/ ead (x,p) =- 



- soft 



2tt 2 Jo 



x 




(B.7) 



But, through the above line and with view to ( C.4 ), the transverse function L^ oft vanishes. 



This completes the derivation of (5D[). 

For the gutter problem of Sec. IV D we must still learn about the first non-vanishing 
piece of Lf {t . Let us work with X t = 1 and remember m — > A*m at the end. We start 
from the full expression, but separate the cut and pole parts of the spectral densities. In 
particular, p t cut ' lead m eans the second term in (|C.2|) , and r 4 lead = 3m 2 /4p the prefactor of 



the delta functions. There is an exact expression (without index lead) to both. Then, 
three differences may be formed : 



2tt 2 



dpp 2 / dx 



1 



+ nix) 



| 2 (V f -r lead ) 6(x-u t ) 



+ Pi 



cut 



-cut, lead 



Pt 



2 

+ 2r/ cad [S(x-(J t )-S{x-p)]} 



(B.8) 



where u t = oj t {p) is the transverse plasma frequency, to be obtained by solving u 2 = 
p 2 + n t (ujt,p). We now notice that x, p are restricted to soft values by the above first two 
differences, but not by the third one. So, in front of the first two, we may still use the 
T/x approximation. Note that \ + n(x) — T/x = /3x/12 + 0{[3 2 x 2 ). Hence, for (3x ~ g 
this difference is by two ^-powers smaller than T/x ~ 1/g. It might contribute to L only 
at g 5 . Working this way, the sum rule helps again to get rid of p t and f t : 

L-(A, = 1) = 

Am 1 Jo \ oj t 

For convenience, this can be further rewritten by introducing u = u t as the integration 
variable (and by once more replacing n(p) — > T/p in a soft term - this time required for 
consistency) : 



/ dp 1 h /?p n (a;*) - n(p) 

JO \ UOt ' 



(B.9) 



Lr (A* = i) 



471-2 



duj [ 1 — (3um(uj) 



UOLO' 



(B.IO) 



31 



with u' the derivative of u t with respect to p, and p being p(u). The square bracket starts 
as for small its saturation at 1 being never reached because the round bracket 
sets the limit. It starts with 1/6 (at uo — > m) and goes as (9/4)m 4 u; _4 ln(c<j/m) for large 
uj (with such details taken from Appendix B of [23|). Hence ( P.10|) is indeed of order g A 



in magnitude. Going to At 7^ 1 simply amounts to m — > Atm in ( |B.10| ). But note that 
this scaling also changes the definition of e.g. uo t , which now is the transverse plasma 
frequency as if m were Atm. 



What we really need in the main text, is not L t itself but the first factor in flBSp. The 
operation there, fortunately, eliminates the above last integration : 



2 L d L - ^ 
Y t Lt ~ dxtLt ~ 216^ 




fB.ll} 



This "first factor" is thus positive, and it behaves as ~ 5 ,4 A^ for small g. Just these 
properties were used in the main text below (|65| ) to reach the long bath-tub. 



Appendix C 



Here we collect a few special details on the spectral densities p t and pi which were 
needed in Appendix B. There we had to learn on the product 



p(x,p) = (x 2 -p 2 ) p(x,p) 
and its asymptotic forms at large p-argument (p 2 3> m 2 ) : 



Pt 



lead 



3m 2 



lead 



4p 
3m 2 



[S(x — p) — 5(x + p) ] 



3m 
4p3 



x 9 (p 



x 



+ 



2p 3 



x 6 (p 



x~ 



(C.l) 

(C.2) 
(C.3) 



These leading terms are readily obtained from the full expressions as given in Appendix 
B of p3[ . One may check ( |C2| ) and ( p.3|) by using it in the p sum rules and thereby 
producing, in each case, the term of highest p -power to the right. The exact p sum rules 
read : 

(1/x) ( 



dx < 



x 



X 



X 



> Ptfap) 



m 

|p < m + m 
^p 4 m 2 + ^■p 2 m A + m 6 

|p6 m 2 _|_ Z2|p4 m 4 _|_ ^p2 m 6 _|_ m 8 



(C.4) 
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( l/x) 



dx < 



x 

3 



X 



X 



X 



> Pt(x,p) 



3m 2 / (3m 2 + p 2 ) 

2 

m 

prm + m 
|p 4 m 2 + |p 2 m 4 + m 6 
\p % m 2 + yi|p 4 m 4 + \p 2 m° + m 8 



3^ 1 175^ 1 5 

They derive through ( |C.1| ) from the sum rules of ordinary densities : 



fl/xl 



dx < 



x 



x 



x 



X 



X 



> Pt{x,P) 



l/p 2 

1 

p + m 
p 4 + yp 2 m 2 + m 4 
£> 6 + ¥rp 4 m 2 + ^-p 2 m 4 + m 6 



x 



X 



X 



X 



X" 



Pt{x,p) 



35^ 1 5 

1297 4™4 
UtsP m 1 5 

l/(3m 2 +p 2 ) 



p 8 + fj|p 6 m 2 + ^p 4 m A + f p 2 m 6 + m c 



p 2 + m 2 
P + %p m + m 



p D + ||p 4 m z + ^p^vri 1 + m D 
p 8 + f§p 6 m 2 + ff|p 4 m 4 + ^p 2 m 6 + m 8 



(C.5) 



(C.6) 



(C.7) 



and these, in turn, are derived along the lines given in f20| 
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